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We study the existence and stability of cosmological scaling solutions of a non-minimally coupled 
scalar field evolving in either an exponential or inverse power law potential. We show that for 
inverse power law potentials there exist scaling solutions the stability of which does not depend on 
the coupling constant f. We then study the more involved case of exponential potentials and show 
that the scalar field will assymptotically behaves as a barotropic fluid when £ <C 1. The general case 

ON ' § ^ 1 is then discussed an we illustrate these results by some numerical examples. 
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I. INTRODUCTION 

On one hand recent cosmological observations, and particuarly the Hubble diagram for Type la supernovae El have 
y—i . led to the idea that the universe may be dominated by a component with negative pressure [0 and thus that today 
^ ' the universe is accelerating. Such a component can also, if one sticks to the prediction of inflation that Q = 1, account 
for the "missing energy". Yet, many candidates have been proposed like a cosmological constant, a "dynamical" 
cosmological constant [|| , cosmic strings [f| or a spatially homogeneous scalar field rolling down a potential || . 

On the other hand, potentials decreasing to zero for infinite value of the field have been shown to appear in particle 
physics models (see e.g. |6|Q]). For instance, exponential potentials arise in high order gravity ||, in Kaluza-Klein 
theories which are compactified to produce the four dimensional observed universe j^] or can arise due to non- 
^\ . perturbative effects such as gaugino condensation fic[ | . Inverse power law potentials can be obtained in models where 
""^ ' supersymmetry is broken through fermion condensates Q|. This gives one more theoretical motivation to study the 
q-( cosmological implications of a field with such potentials. 

The cosmological solutions with such a field have first been studied by Ratra and Peebles [|ll] (see also [ l^,|l3| ) who 
showed the existence and stability of scaling solutions in respectively a field, radiation or matter dominated dominated 
• • ■ universe for a field evolving in an exponential and inverse power law potential. A complete study in the framework 



. £h ', °f barotropic cosmologies in the case of the exponential potential [g_4 1 5 1 show that the solutions were stable to shear 
perturbations and to curvature perturbations when P/p < —1/3, but that for realistic matter (such as dust) these 
solutions were unstable, essentialy to curvature perturbations. Liddle and Scherrer made a complete classification 
of the field potentials and show that power law potentials also lead to scaling solutions (i.e. to solutions such that the 
field energy density p^, behaves as the scale factor at a given power) and the coupling of the field to ordinary matter 
has been considered in [jlTj . Such solutions are indeed of interest in cosmology since they provide a candidate for a 
component with negative pressure. 

Cosmological models with a scalar field have started to be investigated (lqjly] for different kind of potentials as 



cosine potential M, exponential potential |20 2M and inverse power law potentials p2j. It has also been shown that 



the luminosity distance as a function of redshift [|23 24 or the behaviour of density perturbations in the dust era as 
a function of redshift can be used to reconstruct the field potential. 

However, all these studies have been done under the hypothesis that the field is minimally coupled to the metric. 
It is known that terms with such a non-minimal coupling TZf(4>) between the scalar curvature TZ and the field 



can appear when quantizing fields in curved spacetime [E5 26 1, in multi-dimensional theories |27|] like superstring and 



induced gravity theories P8| . Since these theories predict both the existence of scalar fields with potential or power- 
law potential and non-minimal coupling, it is of interest to study the influence of this coupling, and for instance the 
robustness of the existence and stability of scaling solutions. The influence of such a coupling during an inflationary 
period and the existence of inflationnary attractors have yet been examined (see e.g. p9|). 

In this article, we study the stability of scaling solutions of a non-minimally coupled scalar field. We first present 
(§|h|) the main notations and equations. After having, in §111, briefly recalled the standard approach for determining 



the potentials that can give rise to such behaviour for a minimally coupled scalar field, we investigate the inverse 
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power law potentials (§IV) and the exponential potential (§V). In the latter case, we study the two limiting situations 
^ < 1 and £ 3> 1 and then have a heuristic discussion in the general case. 



II. DESCRIPTION OF THE MODEL 



We assume that the universe is described by a Friedmann-Lemaitre model with Euclidean spatial sections so that 
the line element reads 

ds 2 = -dt 2 + a 2 {t) [5 l3 dx l dx j ] = g^dx^dx", (1) 

where a(t) is the scale factor and t the cosmic time. Greek indices are running from to 3 and latin indices from 1 
to 3. The Hubble parameter H is defined as H = a/a, where a dot denotes a derivative with respect to t. We also 
introduce V M the covariant derivative associated to g^ v . 

We assume that the matter content of this universe is composed of a perfect fluid and a homogeneous scalar field 
</> coupled to gravity and coupled to matter only through gravity. The fluid energy density, p B , and pressure, P B , are 
related through the equation of state 

P B =^ B p B , (2) 
where B refers to "background" . The conservation of the energy-momentum of the fluid reduces to 

p B +3H(p B +P B ) = 0. (3) 
The scalar field <fi evolves in a potential V(cj>) and its dynamics is given by the Lagragian 

^ = ~ J %m + 2^Hf(<p) + 2V{<j>)] V^d 4 x, (4) 

where £ is the field-metric coupling constant (£ = corresponds to a minimally coupled field and £ = 1/6 to a 
conformally coupled field), 7Z is the scalar curvature of the spacctime. No known fundamental principle predicts 
the functional form f{4>) and we will assume that f(<p) = <fi 2 /2. This is however the only choice that allows for a 
dimensionless £. The equation of evolution is then obtained by varying the action (^) with respect to the field which 
leads to the Klein-Gordon equation 

^=0~D*-W-a^=0, (5) 

where □ = V^V^. For an homogeneous field in the spacetime (|lj), it reduces to 

+ 3i?0 + ^ + 6£ (W 2 + Hj <j) = 0. (6) 
The energy density and the pressure of this field are defined by 

P<f> = a- 2 T Q0 ; P = (7) 

and we define uj^ = P^j p^,. The energy-momentum tensor is obtained by varying the action with respect to 
the metric g^ v and reads 

Tjt„ = (1 - 2£)V P 0V„0 + ^ - ^.Va^vV - 2e0V^V„0 + 2£0n# M „ + £G^ 2 - V{<j>)g» v , (8) 

being the Einstein tensor of the background metric. This expression can be compared to standard results (see 
e.g. |25f] ). Note however that when the signature of the metric is (+, — , — , — ), one has to change the sign of g^ u , □, 
R^v, and while T^ v and 1Z remain unaffected. 

It is then straigthforward to check that the density and the pressure (0) of the scalar field are given by 

P4> = ltf + V{ct>)+ 3H^(2^+H<f>) (9) 

P4> = \-4> 2 - V{4>) - £ ((2H + 3H 2 )c/) 2 + AH<t>4> + 2^ + 20 2 ) , (10) 
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and that the conservation of the energy-momentum of the field (V^T^ = 0) reduces to the Klein-Gordon equation 
P(t> + 3H( P<f> + P 4> ) = 0^4> + 3Hj)+^- + 6Z (2H 2 +H^4> = 0. (11) 



00 



[We have used that the scalar curvature and the Einstein tensor are respectively given by 11 = 6 (2H 2 + H\G, 
3a 2 H 2 and G l j = —(2H + 3H 2 )a 2 8 l j]. The equation of state of the field is defined by 

^ = ~- (12) 

P<i> 

The matter content being described, we can write the Einstein equations which dictate the dynamics of the spacetime 
and, in our case, reduces to the Friedmann equations 

H 2 = ^(p B +p^), (13) 
H = ~((u B + l)p B +p^ + P <t> ) (14) 

with k = 8itG. One equation of the set (p|jLl], p^JTil ) is redundant due to the Bianchi identities. It will also be usefull 
to introduce the density parameter of a component X as fix = npx /3H 2 . They must satisfy (from equation |l3|) the 
constraint 

fi s +fi* = l. (15) 



III. SCALING SOLUTIONS FOR MINIMALLY COUPLED SCALAR FIELDS (£ = 0) 

In this section, we briefly recall the "standard" procedure to show that there exist scaling solutions and to determine 
the potentials that give rise to such solutions. This presentation will also enable us to understand the differences with 
the more general case of non-minimally coupled scalar fields. We follow the approach initiated by Ratra and Peebles 
jllj and others [ fU3|fl6| ] where one, assuming a scaling form for p^, derives an equation for <p(a) and then uses (p[]lO|) 
to deduce the associated potential. 

The equation (^) of evolution of the background fluid density can be integrated to give 

Pb=Pbo x ~" 1 with x=— , (16) 

do 

where a subscript refers to quantities evaluated at a given initial time. We now look for scaling solutions, i.e. 
solutions such that 

P0 = p<f,QX-' n <S=^ P = Q - l) P4> . (17) 

Since n/3 — 1 = 1 — 2V/ p^ e [—1; 1], we deduce that n £ [0, 6] (let us emphasize that this is a priori no longer true 
when £ 7^ 0). Using equations (p|-|l0j) , such a solution must satisfy 



and V{<j>) = (l- r ^)p <t> . (18) 



n 

-p v . vr , ^ g 



Now, the Friedmann equation (|T3| ) implies that the field should satisfy 



d4> - 1 with B = ,[e™ and A E ^/A (19) 



dx aVl + B 2 x n - m ' V Pd>o V « V P<po + Po V k 

This can be integrated for different relative values of m and n. 
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A. m — n case 



In that situation, (OjJ) leads to 



and then using (fill), to the potential 



^\nx with A^s-^L (20) 

A J UK 
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V^) = [l-'-)p^e- x ^ M . (21) 

This solution corresponds to the scalar field dominated universe of Ratra and Peebles and to their scaling 
solution in the case m = 3 and m — 4 (i.e. radiation or matter dominated universe). Note that with such a potential, 
p§ will, by construction, mimic the evolution of the background fluid and that we do not have to assume that f2^ <C 1. 

Note however that, if the scalar field has reached the attractor from very early time, behaves like radiation and 
thus contributes to a non-negligible part of the radiation content during the nucleosynthesis and it has been shown 
that it implies the constraint f^ < 0.15 20|l^ |. Moreover, since lo^ = in the matter era, such a field will not 



explain the supernovae measurements (which seem to favor uj^ = —0.6 f3C[| ) even if it can account for a substantial 
part of the dark matter. 



B. m 7^ n case 



In this case we have 



A 



dx xVl + B 2 x n ~ 



which can be integrated (B > 0) to give (31j 

, , 2A 



m — n 



■In 



1 + B~ x x- 



B^x^ 



Again, using (|18[), we can deduce the potential 

V(<t>) = (l-^)pwx- n , 



(22) 



(23) 



(24) 



x being given by (|23|). Indeed, we only get the potential in a parametric form, but when B 1 (i.e. when the perfect 
fluid drives the evolution of the universe, x can be eliminated from (p3||24|) to give 



V(4>) = (!-£) P4>v 



2A 



-B 



(<P-<f>oT 



(25) 



When m — 3 and m = 4, we recover the Ratra-Peebles result |TT| as well as the Liddle-Scherrer result jl(| for all m. 
This parametric general form of the potential seems however not to have been exhibited before. 



IV. NON-MINIMALLY COUPLED SCALAR FIELDS WITH A POWER LAW POTENTIAL 



A. Existence of a scaling solution 

The former procedure cannot be applied when the field is non-minimally coupled since it is impossible for instance 
to write a closed equation for ^ as in (jl9|) . Moreover, we are interested in a field evolving in a given potential. 
We assume that the field evolves in an inverse power law potential and show that there exist scaling solutions, the 
stability of which is then studied. 

We assume that the potential takes the form 
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V{<f>) = V M* 



M„ 



with a > 0, 



(26) 



with M p being the Planck mass. The universe is dominated by the perfect fluid so that (we assume m^O) 

mt — iq 3 
Redefining M p (t — to)v / Vo' as t and 4>/M p as cf>, the Klein-Gordon equation (||) becomes 

rat ra \m J * 
Looking for a solution of this equations of the form <j) oc t", one obtains 



(27) 



(28) 



i ,8 ±a+2 
= <Pot ' , <P Q = 



:, = 



so that oc a~ n with n/m = a/(a + 2). This solution is only well-defined if 

6 a 6/4 



m a + 2 ra \ra 



1 (a + 2)f > 0. 



(29) 



(30) 



One can then compute the energy and the pressure of this field by inserting this solution into (p[]l0|) and verify, 
after some algebra, that 



lo b ch — 2 
a + 2 ' 



(31) 



whatever the value of £. This shows that the scaling solution does not depend on the coupling in the sense that ui^ 
is independent of £. This relation generalises the one found for minimally-coupled scalar fields P, pd| , ^6[ . 

B. stability 

As emphasized in the previous section, the scaling solution does not depend on the coupling £. The stability of 
such a solution is known when £ = |ll]]l6| , we now have to study it when £ + 1 0. Following |il|| , we define the new 
set of variables 



t = c T 



u(t) = 



Mr)' 



(32) 
(33) 



where <j> s is the scaling solution (|29|). We set (/)' = ^2. Using <p = c 
equation (pSh reduces to 



-It (All _ ±1 



4f) and <j)'J^s = 2/(a + 2), 



6 



ra a + 2 



-l U' 



a + 2 



6 



to a + 2 



12 

m 



4 

1 

ra 



~(a+l)\ _ 



= 0. 



(34) 



The scaling solution corresponds to the critical point u = 1. Introducing v — u' and linearising around this critical 
point [it = 1 + e] , we obtain 







m a+2 

The eigenvalues, A±, of this system are 



m L m 



± _ 1 £ 1 _ 8. _ _JL_ 

m I ^> I I m a + 2 



(35) 



2A± = 1 



6 



to a + 2 



±W l 



6 



to a + 2 



6 a \ 6 , „, / 4 

+ — a + 2 1 £ 

to a + 2/ to V to 



(36) 
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This expression reduces to the Liddle-Scherrer result |ll| when £ = [with a/ (a + 2) = n/m] and to the Ratra-Peebles 
result 11 when either to = 3 or m — 4. 

A necessary and sufficient condition for the critical point to be stable, is the negativity of the real part of the two 
eigenvalues. Defining £ as 



- (^-6-^) 2 -8m(6-3 
48(a + 2)(4-m) 

the two solutions of (Bq) are real when (£ — D(m — 4) > 0. Thus, 



(37) 



• when the two eigenvalues are complex, A = A^ and their real part is given by, Re(X±) = — ^ — J /2 and 
the scaling solution is stable if 



1-1— i^<0. 
to a + 2 



(38) 



when the two eigenvalues A± are real, their product is (from 



A_A+ = 2 



to a + 2 



6(a + 2) 



1-1 

m 



e >o, 



because of the condition (p0|). They are thus of the same sign, their sum being 

6 4 



A+ + A_ = 1 



to a + 2 

they will both be negative if A + + A_ < and thus, the solution will be stable only if 



G 



to a + 2 



<0, 



as in the case £ > £. 



(39) 



(40) 



(41) 



In conclusion, we find that whatever the coupling constant £, the scaling solution (29) will be stable if and only if 
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1-1 



to a + 2 



< 0. 



(42) 



The value of £ only determines the nature of the stable point, i.e. wether it is a stable spiral or a stable node. This 
generalises the study by Liddle-Scherrer Jll| to a non-minimally coupled scalar field. 

V. NON-MINIMALLY COUPLED SCALAR FIELD IN AN EXPONENTIAL POTENTIAL 



A. Scaling solutions ? 

We now focuse on potentials of the form 

V((f>) = VoM£e- Xct > /M » with A > 



(43) 



and work under the same assumptions as in §1V. Redefining t and cf> as in §[V, the Klein-Gordon equation now reads 

$ + ~-A + - ( - - A tit - Ae"* = 0. (44) 

TO I TO y TO / t 

When the coeficient of cf> vanishes, i.e. in a radiation domi nated universe (to = 4) or when £ = 0, one can find a 



special solution of the form <f> = In (At"). We get that (see § [II A with n — M p 2>1 



G 



s =ln(A^); A- A = A^_iV £ = ^ = ™ _ !, (45) 

where the subscript s refers to the scaling solution. In radiation era, this implies for instance that uj ( p g = 1/3 and the 
scalar field behaves like radiation. Indeed, this is a very special case since in that period TZ — and the field does not 
"feel" the non-minimal coupling and it evolves as if it were minimally coupled. The complete study of these solution 
in function of the two parameters (A, to) @,^5| shows that when A 2 > to, the scaling solution (j> s is a stable node or 
spiral whereas, when A 2 < to, the late time attractor is the field dominated solution, which we do not consider here. 
The convergence towards the solution <p s is illustrated on figure ^. 

Now, in the most general case where to 7^ 4, it is easy to realize that a solution of the form cj> = In (At 13 ) cannot be 
solution of equation (^) . 



1. The |£| < 1 case 



Let us first look at the effect of a small perturbation in £ in the sense that the potential term dominates over the 
coupling term in the Klein-Gordon equation. For that purpose, we set 



u = In t 

^ = ^ + ^ + o(e 2 )- 

The equation of evolution for tp can be deduced from the Klein- Gordon equation (p4|) 



G 



rn 



12 



i/>" + - - 1 W>' + - - - 1 = -- [1 



TO \ TO 



A 



TO \ TO 



(46) 



(47) 



where a prime denotes a derivative with respect to u. Now, if we restrict to £ <C 1 and linearise this equation using 
the expression ( fi5| ) for cf> s , we obtain at zeroth order in £ 



i-lH^lA-l),."^-!)^ 



4 

TO 



-A 



2u], 



the solution of which has the general form 

ifj = B+e a+u + B^e a ~ u 



12 4 - to 
toA 6 — to 



1 + lnA 



—A 



(48) 



(49) 



B + e a+u and B_e a - U are two independent solutions of the homogeneous equation. Since a + a_ = 2(6/m — 1) = 
— 2(a_|_ + £*-)) we deduce that if to < 6, the real parts of both a+ and a_ are negative, so that the two homogeneous 
solutions correspond to decaying modes and the particular solution is then an attractor (if to < 2/3 then a + and a_ 
are real, otherwise they are complex and the homogeneous part will decay while oscillating). 

Indeed, this analysis is valid only as long as the potential term dominates over the coupling term in equation (|44|), 
that is, as long as 



dV 



> 



TO \ TO 



-2d 



u <C u e 



_ 1 to(6 — to) 1 



12|f| \A-m\ A 



In- W2 1 



G 



(50) 



dY 



00 and we are back to the minimally coupled case (§111 A). The 
and ^ (^ — l) £0 in function of u is shown on figure ||. When £ -C 1, we see that, as expected, the 



We recover that when £— > or m ^ 4, 11 
behaviour of 

solution is first dominated by the potential term but that, as time elapses the coupling term tends to become more 
dominant. In figure we illustrate how the phase space trajectories are deformed due to the existence of this small 
coupling. 

Now, as long asu< u eq , we can compute the equation of state of the field by inserting the particular part of the 
solution (E9|) in (&[L0|) which will take the form 



W0(ti) = u^ s (1 + A(to,A)u" 1 +S(to,A,£)m" 2 + 0(u^ 3 ) if to^3 
= --(1+C(to,A)u" 1 +C'(u- 2 )) if to = 3. 



(51) 
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The difference in these two behaviours comes from the fact that P^ oc u 2 if to ^ 3 and P^ oc u when m = 3. The exact 
forms of the functions A, B and C can be obtained by doing an expansion of P<p/ P4, in u^ 1 . In figure [| we show the 
deviation of the equation of state from pure scaling and compare the former expansion to the numerical integration. 
In conclusion we have that when ( < 1 and u -C u eq , the field converges towards a barotropic fluid. 



2. The |£| > 1 case 

Let us now consider the case where initially the coupling term dominates over the potential term in the Klein-Gordon 
equation. At lowest order, one has 

i-lV + -(±-lU^0, (52) 
.to ) m V to 

the solution of which is of the form 

4> = A+c a + u + A_e a - U . (53) 

When £ > 0, since a + a- oc 4/m — 1 and a+ + a_ = — (6/m — 1), we deduce that if < to < 4, the real part of both 
a + and a_ are negative so that (53) corresponds to two decaying modes [It can be seen as a proof that the critical 
point <p = is an attractor] so that 

ll^ct) -> and Ae~ A0 -> A when < m < 4 and £ > 0, (54) 

and the potential term will rapidly catch up the coupling term. 

When £ < 0, a + a- oc — 4/m+ 1 and the real part of one of the two quantities, a+ say, is positive when < m < 4 
so that 

7^0 -> ex) and Ae~ A0 -> when < to < 4 and £ < 0, (55) 
and the coupling term will dominates forever (see figure 0) and cf> will behave as A + e a+u with 



^ = i.-i) + , ^-iV-^(i-i]f, 

TO/ V\Ti / TO\TO ' 



for which, since the potential term is negligible, 



2 



4£ [£(3-m) + (£-l) «+ + 2 C 



This solution has however to be excluded since one can check that it leads to p& < 0. 



3. General case 

The general case is more involved since we cannot find any analytic solution to (p4|). First, when £ > 0, we have 
seen that when either the potential or the coupling term dominates, the other slowly catch up. We thus expect a late 
time solution which satisfies 

K£<j> ~ Ae" A0 . (58) 

In figure |], we plot the evolution of these two terms in the case where ( < 1 and in a more general case. The two 
terms alternatively dominate and then converge to the same value. Indeed, this is no proof. 

We can however look for a general solution of the form 

n— 00 
n=0 
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with tfjQ given by ( flq ) and if>i given by ([49]). Inserting this expansion in the Klein-Gordon equation (|44|); we obtain 
the hierarchy of equations 

K +(--l) V4 + 2(--l) V>„ = -—(-- l) ^ n -i + AA- A /n(^o, .... Vn-i) if n > 1. (60) 



m \m 



The functions /„ depends on all the solutions ^ for i < n. All these equations have a solution of the form 

if) = B+c a + u + B_c a - U + i> n . (61) 



As in §VA2, the two homogeneous modes decay if m < 6 and Vn,^ — » ^>„. It can also be seen that ipn will be a 
polynomial in u of degre n. If the series £""0" converges then, 

£™t/i„ when u — > oo, (62) 



n— oo 



n=0 



from which we can conclude that if X^C'Vn converges, there exists an attractor to the equation (Q) given by (|62j). 
Indeed, we cannot demonstrate the convergence of this series in the general case. In figure || we show the phase space 
trajectories showing the convergence towards this attractor. 

When m > 2/3, the solution converges towards the attractor while oscillating so that the equation of state have 
wiggles and converges towards ui ( f >s (see figure^). When m < 2/3, this oscillatory behaviour does not appear whatever 
£ (see figure ||) . 



4- Case of a radiation dominated universe 



As explained above, one has the solution ([15|) to the Klein-Gordon equation and this solution is an attractor (see 
the phase analysis of figure |5|) and the field behaves as radiation. The solution ( p9| ) reduces to the two decaying 
modes. Indeed the derivation of (|5l]) is no longer valid but, with the same method, one can show that, to first order 
in u~ l , 

^ = i(l+-4(A) U - 1 +6(A,e)^ 2 ), (63) 



where the two functions A and B are obained as in § V A 2 . Note that for u 3> 1 the £ terms of (|9f]l0|) dominates the 
density and the pressure. In figure ^|, we have plotted the phase space showing the convergence towards the attractor 
<fi s in the case of a conformally coupled scalar field and the evolution of its equation of state. 

B. Numerical results 

We integrate numerically equation ( ff4| ) and use the Poincare projection 63] to represent the result in the plane 

(0,0') 

r 

■ sin 6* 



1 - r 

cj)' = — cos 9. (64) 
1 — r 

This projection shrinks all the trajectory in the phase-space to the unit disk. The point N = (0,1), E = (—1,0) 
and W = (1,0) respectively correspond to ((f) = oo, <f>' = 0), (4> = 0, <\>' = — oo) and ((f) = 0, (f)' = oo). Note that the 
system is non-autonomous since equation ( pl7| ) depends on u. Thus, curves can cross in the Poincare representation 
but they will not cross in a 3D representation ((f>, 4>',u). The following plots corresponds to the situation studied in 
the previous paragraphs. We plot the evolution of the solution in the phase space (</>, <f>') and the evolution of the 
equation of state when £ <C 1. To finish, we give some examples of evolution in the case £ ^ 1. 








m=3, h=2, |=0.01 m=3,h=2, £=0.1 m=3, X=2, £=-0.01 




FIG. 2. The deformation of the phase space from figure [l] due to the coupling for £ = 10 2 , 10 \ — 10 2 . The solution 
converges towards the attractor which drifts from cj> s , since <j> — <j> s oc £u + 0(£ 2 ). 



m=3, A=2 m=3j x=2 




FIG. 3. The time varying equation of state and its deviation from the minimally coupling case for £ = 0, 10 3 , 10 2 and the 
comparaison of this equation of state with the estimate ( |5l| ) at first order and second order in u^ 1 . 
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m=3, X=2, £=1/2 



m=3, X=2, £=10"' 



m=3, h=2, (,=-0.01 



■ \ 
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FIG. 4. The respective influence of the coupling term and the potential term in the Klein-Gordon equation. We have plotted 
|^| and |7££0| normalised to their sum for £ = 1/2 (left). The two terms alternatively dominate and then converge. When 
^Cl, after some oscillations corresponding to the convergence towards <j> a , the potential term dominates over the coupling as 
long as ti< u eq . When £ < 0, the coupling term will dominates forever. 





VI. CONCLUSION 



In this article, we studied the influence of the coupling between the scalar curvature and the scalar field on the 
existence and stability of scaling solution of this field evolving in either an exponential or an inverse power law 
potential. The motivation for considering such solutions is first that they can be a candidate for a matter component 
with negative pressure and second that they appear for a large class of potentials predicted by some theories of hight 
energy physics. 



11 



We first found a new parametric form of the potential that reduces to the inverse power law potential when the 
energy density of the fluid drives the evolution of the spacetime. 

Concerning the inverse power law potentials, we showed analytically that the existence and stability of a scaling 
solution does not depend on the coupling £ and that the equation of state of the field was always given by — "„" 2 . 
This generalises the work by Ratra-Peebles and Liddle-Scherrer |ll| . 

The situation is more involved with exponential potentials since one cannot find an analytic form for a scaling 
solution (apart for a radiation dominated universe). We then studied the effect of a small perturbation coupled to the 
scalar curvature and computed lo^, to first order in £ -C 1 and showed that there always exists a time after which one 
cannot neglect the effect of the coupling. In that limit, we show that the equation of state was converging towards a 
barotropic equation of state and some numerical examples tend to show us that it should be the case whatever £ > 
(but this has not been demonstrated). Indeed, such potentials are not the most favored since they are constrained 
by nucleosynthesis to fl^Q < 0.15 ]20| and cannot explain (when £ = 0) the supernovae measurements since uj^ = 0. 
Note however that the convergence toward the barotropic equation uo^^ is much slower when £ ^ 0. When £ < 0, we 
have shown that there always exists a time after which the coupling term will dominate and thus that there exists a 
scaling solution but with a different equation of state (as long as < m < 4, and thus in a matter dominated era). 
Unfortunately, such a solution has to be rejected since it has negative energy. 

As pointed out by Caldwell et al. [^2|, a smooth time dependent field is unphysical since "one has to take into 
account the back-reaction of the fluctuations in the matter components" . The cosmological consequences of such an 
inhomogeneous coupled scalar field in both exponential and inverse power law potentials (such as the computation 
of the CMB anisotropics and the matter power spectrum) will be presented later Q and our present study is only 
related to the implication of the homogeneous part of such a field. 
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